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Using muon spin resonance we examine the organometallic hybrid compound Cu(l,3- 
benzenedicarboxylate) [Cu(l,3-bdc)], which has structurally perfect spin 1/2 copper kagome planes 
separated by pure organic linkers. This compound has antiferromagnetic interactions with Curie- 
Weiss temperature of —33 K. We found slowing down of spin fluctuations starting at T = 1.8 K, and 
that the state at T ^ is quasi-static with no long-range order and extremely slow spin fluctuations 
at a rate of 3.6 ^sec~^. This indicates that Cu(l,3-bdc) behaves as expected from a kagome magnet 
and could serve as a model kagome compound. 



The experimental search for an ideal two dimensional, 
spin 1/2, kagome compound, which has no out-of-plane 
interactions and no impurities on the kagome plane, has 
powered tremendous experimental effort in recent years 
Yet, all compounds studied today have shortcom- 
ings. Recently, a promising copper-based metal organic 
hybrid compound Cu(l,3-bdc) was synthesized by Nytko 
et al. Q. This compound has an ideal kagome lattice 
structure as indicated by X-ray, the spins are naturally 
1/2, and there are no Zn ions or any other candidates 
to substitute the Cu on the kagome plane. The goal of 
this paper is to show that from a magnetic point of view 
Cu(l,3-bdc) shows the signatures of the high degree of 
frustration expected on the kagome lattice. This is done 
by demonstrating that the inter-plane interactions are 
small enough compared to the intra-plane interactions 
that no long range order is found at temperatures well 
below the interaction energy scale, and by characteriz- 
ing the ground state properties. The experimental tool 
is muon spin resonance (/iSR). Our major finding is that 
the state at the lowest temperature investigated is quasi- 
static with extremely slow spin fluctuations. This type of 
behavior is similar to a huge class of frustrated magnets 
0,13. Therefore, Cu(l,3-bdc) could serve as a model spin 
1/2 kagome compound. 

Cu(l,3-bdc) is shorthand for Cu(l,3- 
benzenedicarboxylate). The kagome planes are separated 
by organic linkers, each linker being a benzene molecule 
with two corners featuring a carboxylate ion instead of 
the standard H ion. If one were to label the corners 1-6 
consecutively, the two corners with the carboxylate ions 
would be the 1st and 3rd. The Cu ions located on the 
kagome plane are linked by 0-C-O ions, while inter-plane 
Cu ions are linked by 0-5C-0 ions. The basic elements 
of Cu(l,3-bdc) are depicted in Fig. [1] Magnetization 
measurements found antiferromagnetic Curie- Weiss 
temperature of -33 K Strong antiferromagnetic 



exchange between Cu^+ ions linked by a carboxylate 
molecule was also found in the trinuclear compound 
Cu3(02Ci6H23)-1.2C6Hi2 [5j]. Heat capacity shows a 
peak at T ~ 2 K tl. 

The powder we examined contains Cu(l,3-bdc) in the 
form of blue crystalline plates. However, it is mixed with 
some green plates of copper-containing ligand oxidation 
by-product C32H24CU6O26, which cannot be separated 
from the blue plates. To the naked eye it looks as if about 
10% of the plates are green. However, as we demonstrate 
below, the ^SR signal from Cu(l,3-bdc) can be separated 
from the C32H24CU6O26 signal. Our sample was pressed 
into a Cu holder for good thermal contact. 




FIG. 1: (Color online) The Cu(l,3-bdc) structure showing the 
kagome planes and the inter-plane and intra-plane superex- 
change path. 



Muon spin rotation and relaxation (/iSR) measure- 
ments were performed at the Paul Scherrer Institute, 
Switzerland (PSI) in the low temperature facility spec- 
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trometer with a dilution refrigerator. The measurements 
were carried out with the muon spin tilted at 45° rela- 
tive to the beam direction. Positrons emitted from the 
muon decay were collected simultaneously in the forward- 
backward (longitudinal) and the up-down (transverse) 
detectors with respect to the beam direction. Transverse 
field (TF) measurements, where the field is perpendicular 
to the muon spin direction, were taken at temperatures 
ranging from 0.9 K to 6.0 K with a constant applied 
field of H = 1 kOe. Zero-field (ZF) measurements were 
taken in the longitudinal configuration at a temperature 
ranging from 0.9 K to 2.8 K. The longitudinal-field (LF) 
measurements, where the field is parallel to the muon 
spin direction, were taken at several different fields be- 
tween 50 Oe and 3.2 kOe with a constant temperature of 
0.9 K. We also performed a field calibration measurement 
using a blank silver plate providing the muon rotation 
frequency = 13.67 MHz at the apphcd TF of 1 kOe. 



in both insets of Fig. [2] The broad and disappearing peak 
in the main panel corresponds to the fast decaying signal 
in the first 0.2 /isec seen in the lower inset. The arrow in 
the inset demonstrates the frequency shift. Consequently 
we fit the function 

Arpit) = Aie-^"''*^''/^cos{ujit+(p)+A2e-^^'^cos{uj2t+(p)+Bg 

(1) 

to our data in the time domain globally, where the param- 
eters Ri and wi are the relaxation and angular frequency 
of the by-product, and R2 and uj2 are the relaxation and 
angular frequency of the kagome part. The parameters 
Ai = 0.049(4), A2 = 0.125(3), i?i = 0.13(1) (Atsec'i), 
if and Bg are shared in the fit, while R2tUJi,uj2 are free. 
The quality of the fit is represented in the inset of Fig. [2] 
by the solid lines. The ratio of Ai to A2 supports the 
assignment of the fast relaxing signal to Cu(l,3-bdc). 
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FIG. 2: (Color online) FFT of the asymmetry data in a field 
of 1 kOe transverse to the initial muon spin dirrection. fs is 
the reference fequency in pure silver. Inset: transverse field 
asymmetry in the time domain and a rotating reference frame. 



In the inset of Fig. [2] we depict by symbols the 
muon decay asymmetry in a reference frame rotated at 
H = 200 Oe less than the TF. In the main panel of 
Fig. [21 we show the fast Fourier transform (FFT) of the 
TF data at some selected temperatures. The FFT of the 
highest temperature, 6 K, shows a wide asymmetric peak 
with extra weight towards low frequencies. At 3 K the 
wide asymmetric peak separates into two different peaks 
shifting in opposite directions. At even lower tempera- 
tures the low frequency peak vanishes. We assign the 
latter peak to muons that stop in Cu(l,3-bdc) since such 
a wipe-out of the signal is typical of slowing down of spin 
fluctuations, which in Cu(l,3-bdc) is expected near 2 K. 

Despite the disappearance of the second peak in the 
frequency domain, its contribution in the time domain 
is clear. The high frequency peak in the main panel of 
Fig. [5] corresponds to the signal surviving for a long time 




FIG. 3: (Color online) The muon shift from C32H24CU6O26 
Ki and from Cu(l,3-bdc) K2, and the relaxation rate from 
Cu(l,3-bdc) R2 versus temperature. 



In Fig.[3]we plot the shift, Ki^2 = (wg— wi,2)/ws, versus 
temperature, where tOs = 27r/s. As expected K2 increases 
with decreasing temperatures. The small decrease of Ki 
is not expected and is not clear to us at the moment. 
The muon transverse relaxation, R2, is also presented in 
Fig. [51 It has roughly the same temperature behavior as 
the shift, K2- However, at T =1.8 K R2 seems to flatten 
out before increasing again around 1 K. This is somewhat 
surprising. 

The field at the muon site is given by B = H — 
J2k (Sfc) where (Sfc) is the thermal average of the spins 
neighboring the muon, and Ak is the hyperfine interac- 
tion with each neighboring spin. Assuming a distribution 
of hyperfine fields in the z direction one can write Ak as 
a sum of a mean value Ak plus a fluctuating component 
SAk- For the distribution 



p{SAk) = - 



CTk 



1 



^{6Akf+al 



(2) 
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one finds that the shift is riven by K2 = (Sfc) J^k ^k/H 
and i?2 = Ifj. (S/c) o'fc is]- To test this derivation it 
is customary to plot the macroscopic magnetization M 
measured with a SQUID magnetometer versus i^2. This 
is depicted in Fig. IJlJa). The magnetization is also mea- 
sured at 1 kOe. The plot indicates that in the temper- 
ature range where both M and K2 are available they 
are proportional to each other. Therefore, Ak is 
temperature-independent . 

If the ak are also temperature- independent parameters 
we expect R2 ex. A plot of R2 versus shown in 
Fig-IHb), indicates that R2 is not proportional to or even 
does not depend linearly on K2 and a kink is observed at 
To = 1.8 K. This result suggests a change in the hyperfine 
fields distribution at Tq. An interesting possible explana- 
tion for such a change is a response of the lattice to the 
magnetic interactions via a magnetoelastic coupling Q. 
However, unlike a similar situation in a pyrochlore lattice 
Q, it seems that here the lattice is becoming more or- 
dered upon cooling since the rate of growth of R2 below 
Tq is lower than at higher temperatures. 
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FIG. 4: (Color online) (a) Magnetization versus the muon 
shift in Cu(l,3-bdc). (b) The muon relaxation rate versus the 
muon shift in Cu(l,3-bdc). 



tion. However, the field distribution is so wide that the 
oscillation is damped quickly. The origin of the recov- 
ery is the fact that some of the muons experience nearly 
static field in their initial field direction during the entire 
measured time. These muons do not lose their polariza- 
tion while others do. When the external field increases, 
the dip moves to earlier times (as the field scale increases) 
and the asymptotic value of the asymmetry increases as 
well (as more muons do not relax). 

The ZF data at three different temperature are shown 
in Fig.[5l^b). As the temperature decreases the relaxation 
rate increases due to the slowing down of spin fiuctua- 
tions, until at the lowest temperature the dip appears. 
We saw no difference in the raw data between 1.0 and 
0.9 K and therefore did not cool any further. 

These are unusual /iSR data in a kagome magnet, in 
the sense that the spin fluctuations are slow enough com- 
pared to the internal field scale to expose the static na- 
ture of the muon spin relaxation function, namely, the 
dip, and to allow calibration of the internal field distri- 
bution. Other kagome magnets show the same general 
behavior but without this dip 0]. The data indicate the 
absence of long-range order and the presence of quasi- 
static field fluctuation. If the ground state had long-range 
order, the muon would have oscillated several times due 
to the internal magnetic fleld. Similarly, if the ground 
state was dynamic we would not have seen a recovery of 
the muon polarization after a long time. 
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The /iSR LF data including ZF are presented in Fig. [5l 
The LF data at the lowest temperature of 0.9 K are de- 
picted in panel (a). At this temperature and a field of 50 
Oe, the muon asymmetry shows a minimum at around 
0.1 /xsec. At longer times the asymmetry recovers. The 
origin of this dip is the presence of a typical field scale 
around which the muon spin nearly completes an oscilla- 



FIG. 5; (Color online) The asymmetry at various longitudinal 
fields and T=0.9 K (a). The asymmetry at zero field and 
various temperatures (b) 

To analyze this type of muon spin relaxation func- 
tion, a theoretical polarization function P{h',A,H,t) 
must be generated. It depends on the random field 
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distribution p(B), the spin fluctuations rate v defined 
by (B(i)B(O)) = (S^) e'"^"* , where B is the internal 
local field, and the LF H. In ZF or small LF field, 
standard perturbation methods for calculating relaxation 
functions do not apply and a special method for calculat- 
ing P is required. This function is produced in two steps. 
In the first step the static muon polarization is generated 
using the double projection expression 



P,(0,A,F,t)= / p(B 



|| + ^^cos(7jB|t) 



We found that the Gaussian field distribution 



p(B) = 



(27r)3/2A3 



exp 



2A2 



d^B 

(3) 

(4) 



works best. This Pz{Q,/^,H^t) is known as the Static 
Gaussian Kubo-Toyabe LF relaxation function p^ . 

In the second step the dynamic fluctuations are in- 
troduced. One method of doing so is using the Voltera 
equation of the second kind [11] 

V / dt'P^{i^,H,A,t-t')e-''*'P^{0,H,A,t') 
Jo 

(5) 

The function Pz{0, H, A,t') is taken from the first step. 
The factor e^''* is the probability to have no field changes 
up to time t. The factor e"''* vdt' is the probability den- 
sity to experience a field change only between t' and 
t' -t- dt' . The first term on the r.h.s is the polarization 
at time t due to muons that did not experience any field 
changes. The second term on the r.h.s is the contribu- 
tion from those muons that experienced their first field 
change at time t' . The factor e-^^'P^O, H, A,t')iydt' is 
the amplitude for the polarization function evolving from 
time t' to t, which can include more field changes recur- 
sively. This equation can be solved numerically 9] and 
Pz{v^ H, A, t) is known as the Dynamic Gaussian Kubo- 
Toyabe LF relaxation function [10]. 

The experimental asymmetry is fitted with A^p = 
AqPz{v,H, A,t) + Bg. The relaxation from the second 
green phase is very small and is absorbed in the back- 
ground factor Bg. In the fit of the field-dependence 
experiment at the lowest temperature, presented in 
Fig. El^a) by the solid lines, A,i',Aa and Bg are shared 
parameters. We found A — 19.8(4) MHz and v — 
3.6(2) /isec~^. This indicates that the spins are not com- 
pletely frozen even at the lowest temperature. 

When analyzing the ZF data at a variety of temper- 
atures, shown in Fig. [5l[b) by the solid lines, we permit 
only V to vary. The fit is good at the low temperatures 
but does not capture the 2.8 K data at early times ac- 
curately. However, the discrepancy is not big enough to 
justify adding more fit parameters. We plot the tem- 
perature dependence of the fiuctuation rate in Fig. [S] 



f hardly changes while the temperature decreases from 
T 2.8 K down to Tq = 1.8 K. From Tq, decreases with 
decreasing temperatures, but saturates below 1 K. This 
type of behavior was observed in a variety of frustrated 
kagome (Ref. i) and pyrochlore (Ref. |^]) lattices. It is 
somewhat different from classical numerical simulations 
where v decreases with no saturation [l^, E^l- In f&ct, 
the numerical v is a linear function of the temperature 
over three orders of magnitude in T [l^ . 

The inset of Fig.[5]shows v as a function of temperature 
' near Tq on a log-log scale where slowing down begins. 
Only near Tq are our data consistent with a linear relation 

V -Voo = vq{T - To), 

where v^o is the high temperature fiuctuation rate. The 
discrepancy with the numerical work might be because 
/iSR probes field correlations involving several spins near- 
ing the muon, while the simulations concentrate on spin- 
spin auto correlations (with a decay Ta compared here 
with v). At our lowest temperature the rotations of en- 
semble of spins are already coherent therefore field and 
spin correlations are not identical. Another possibility is 
that the saturation of v with decreasing T is a pure quan- 
tum effect not captured by the classical simulations. 
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FIG. 6: (Color online) The fluctuation rate u versus temper- 
ature. Inset; v near 1.8 K on a log-log scale. The error bars 
are smaller than the symbol size. 



To summarize, we found that Cu(l,3-bdc) has a special 
temperature Tq = 1.8 K. Upon cooling, the susceptibil- 
ity, as measured by the /iSR, grows monotonically even 
past this temperature. The muon spin line-width also 
grows but halts around this temperature. This might be 
explained by a subtle structural transition, but low tem- 
perature structural data are required. At Tq the slowing 
down of spin fluctuations begins, but the spins remain 
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dynamic with no long range order. The rate of the spin 
fluctuations appears to be hnear near Tq, but becomes 
saturated at the lowest T. This general behavior is simi- 
lar to other kagome compounds, though new features are 
seen here. Therefore, considering its lattice, Cu(l,3-bdc) 
could serve as a model compound for spin 1/2 kagome 
magnet. 
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